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Abstract 

The transverse response function RT{q,u!) for ^He is calculated using the configuration space 
BonnA nucleon-nucleon potential, the Tucson-Melbourne three-body force, and the Coulomb po- 
tential. Final states arc completely taken into account via the Lorentz integral transform technique. 
Non-relativistic one-body currents plus two-body tt- and p-meson exchange currents as well as the 
Siegert operator are included. The response Rt is calculated for q=17A, 250, 400, and 500 MeV/c 
and in the threshold region at q=17A, 324, and 487 MeV/c. Strong MEC effects are found in low- 
and high-energy tails, but due to MEC there are also moderate enhancements of the quasi-elastic 
peak (6%-10%). The calculation is performed both directly and via transformation of electric 
multipoles to a form that involves the charge operator. The contribution of the latter operator is 
suppressed in and below the quasielastic peak while at higher energies the charge operator repre- 
sents almost the whole MEC contribution at the lowest q value. The effect of the Coulomb force 
in the final state interaction is investigated for the threshold region at g=174 MeV/c. Its neglect 
enhances Rt by more than 10% in the range up to 2 MeV above threshold. In comparison to 
experimental data one finds relatively good agreement at q=250 and 400 MeV/c, while at q=500 
MeV/c, presumably due to relativistic effects, the theoretical quasi-elastic peak position is shifted 
to somewhat higher energies. The strong MEC contributions in the threshold region are nicely 
confirmed by data at q=324: and 487 MeV/c. 

PACS numbers: 25.30.Fj, 21.45.-v, 21.30.-x 
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I. INTRODUCTION 



Electromagnetic interactions in the trinucleon systems play an important role in testing 
NN and 3N forces as well as nucleonic current operators. Among the many reaction observ- 
ables available are the response functions which determine the inclusive electron scattering 
cross-section 
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where uj is the electron energy loss, q is the magnitude of the electron momentum transfer, 
is the electron scattering angle and Q = {q, = q^ — '^'^ ■ RL{q,^^) and Rx^q^oj) are 
called the longitudinal and transverse response functions respectively. In order to calculate 
either of these response functions one needs to be able to take into account all final states 
(usually in the continuum) which are connected to the ground state via the current or charge 
operators. This can be accomplished in several ways. For example recent calculations of 
these response functions by Golak et al [l| and Deltuva et al [2] base their calculations on 
Faddeev techniques while we employ the Lorentz integral transform (LIT) j^, 4] method. A 



recent review article on the LIT approach is given in 



The longitudinal response is driven by the nuclear charge density operator and has re- 



cently been calculated covering large parts of the non-relativistic regime in Refs. 
It is notable that although some of these groups use considerably different calculational 
techniques they obtain similar results for Ri{q,uj). All these non-relativistic calculations 
of RL{q,uj) show quite good agreement with experiment for modest momentum transfers 
i.e. q < 400 MeV/c. However for larger q the position of the quasi-elastic peak is sensitive 
to relativistic corrections in the kinetic energy. In the non-relativistic calculation was 
extended up to g=700 MeV/c by choosing a proper reference frame, where relativistic effects 
on the kinetic energy are minimized. A subsequent transformation of the theoretical results 
to the laboratory frame led in fact to a much better agreement with data. Concerning the 
realistic NN interaction model there appears to be a relative insensitivity to which model is 
used. As far as 3N forces are concerned there is no unique picture. For ^He their inclusion 
improves the agreement with data, whereas for one observes the opposite effect 

In RT{q,uj) it is the nuclear transverse current density which drives the response. This 
current density can be expressed as the sum of various components: the normal one-body 
currents with their relativistic corrections, two-body currents arising from meson exchange 



NN forces (MEC) and isobar excitations, three-body currents arising from NNN forces. In 

i n n 

jj the AV18 NN potential [7[ with the UrbanalX NNN potential [8[ were used and one-body 
currents as well as vr- and p-MEC were taken into account to calculate RT{q,uj) at q=200, 
300, 400, and 500 MeV/c and the low-energy Rt at various q. In [2] the CD-Bonn potential 
and its coupled channel extension CD-Bonn+A was taken. The one-body current, vr- 
and p-MEC and A-currents were considered. In addition to computing Rt at q=300 and 
500 MeV/c near threshold responses at various q were also shown. 

Here we present the first fully realistic computation of Rt with the LIT method (in [IC 
the LIT method was applied to the Rt of '^He but with approximations for the MEC and 
using a semirealistic NN potential only). We use the configuration space BonnA potential 



ll| (hereinafter referred to as BonnRA) together with the TM' 12] NNN potential to 



calculate the response at q=174, 250, 400, 500 MeV/c and at q=324 and 487 MeV/c in the 
near threshold region. Our reason for choosing the BonnRA pontential is that the MECs 
are uniquely defined in the case of a boson exchange based potential. We would like to 
emphasize that the LIT method allows us to include consistently the Coulomb interaction 

n n 

in initial and final states, which is not done in UJ, yl. For the electromagnetic current 

nrf 

operator we include the non-relativistic one-body operators plus, as in [1|, 12|], the vr and p 
two-body MEC currents. As is well known the MEC are intimately connected to details 
of the Hamiltonian through the requirement of charge conservation. With boson-exchange 
potentials like the BonnRA and partially for CD-Bonn (contains two effective cx-mesons 
with partial-wave dependent parameters) the form of the MEC are determined by explicit 
knowledge of the boson-nucleon coupling. For phenomenological potentials such as the 
AVIS which was used in [l| one can construct a consistent vr- and p-MEC 15, 16, 17] by 
interpreting the isovector part of the given potential model as due to an effective vr and p 
exchange. 

Our calculation is performed in two ways depending on how we treat the electric mul- 
tipole operators. One method, which we refer to as the direct method, simply uses the 
current operators per se in the electric multipoles. In the second method the electric multi- 
pole operators are transformed via use of the continuity equation into a form which includes 
the charge operator. We refer to this latter form of the electric multipole operator as the 
Siegert form. If both the continuity equation were fulfilled exactly and dynamic equations 
were solved exactly then these two ways would lead to the same results. Since, as in our 
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case, a realistic nuclear force includes components additional to one-boson exchange poten- 
tials, such as momentum-dependent NN forces and 3N forces, the continuity equation is 
only approximately fulfilled when one employs only the dominant, well established MECs. 
Therefore performing calculations in the two ways allows us on one hand to find out to 
what extent the vr and p exchange currents we use are compatible with the realistic nuclear 
force employed. On the other hand, via use of the charge operator it permits us to take 
into account a part of the additonal MEG thus checking their possible relevance. In 
such an investigation has not been carried out (in {l| the Siegert operator is only used for 
reactions with real photons). 



II. NUCLEAR FORCES AND THE CURRENT OPERATOR 

The transverse response Rt which depends on the transverse nuclear current density 
operator Jt is given by 

Here Mt is the mass of the target nucleus, \I^o and denote the ground and final states, 
respectively, while Eq and Ej are their eigenenergies, 

(/i-Eo)^o = 0, {h-Ef)^f = Q, (3) 

where h is the intrinsic nuclear non-relativistic Hamiltonian. States of our system are repre- 
sented by products of normalized center of mass plane waves ipiV^j) and internal substates 
\l/o,/ entering ([2]). Correspondingly, the current operator J in ([21) is related to the primary 
current operator J as follows, 

J5(P;-Po-q) = MP/)|J|y.(Po)), (4) 

where the matrix element is defined in the center of mass subspace. The cross section we 
need corresponds to the laboratory reference frame and we set in (jl]) Pq = 0. The quantity 

is that component of J which is orthogonal to q. The second summation (integration) 
in ([2]) goes over all final states belonging to the same energy Ej, and Mq is the projection 
of the ground state angular momentum Jq. 

The Hamiltonian h includes the kinetic energy terms, the 2N and 3N force terms, and 
the proton Coulomb interaction term. As in [^] the ground state \l/o is calculated via an 
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expansion in basis functions which are correlated sums of products of hyperradial functions, 
hyperspherical harmonics and spin-isospin functions. In the present work the 2N + 3N 
interactions are taken as the Coulomb+ BonnRA+TM' (A=2.835 fm~^) as in ^. The TM' 
cut-off parameter A properly fixes the binding energy to 8.47 MeV. 

We perform a non-relativistic calculation. The current J includes one-body and two-body 
operators. The one-body current operator as obtained from (jl]) is 

A 
k=l 

where A is the number of nucleons in the target nucleus and 

j(fc)spm = e'l GM{k), 

m, = e^""-''' ^ GE{k). 

Here = — Rem, Pa; = Pfc ~ PcmM, and ak are the relative coordinate, momentum, 
and spin operator of the k-th particle and M denotes the nucleon mass, while Rem and 
Pcm are the center of mass coordinate and momentum variables of the A-body system. The 
component jg does not contribute to J^. However, separate multipoles as defined below 
depend on this component. 

In the above expressions we use the notation 

GEMk) = G%^M')^-^ + Gl^M')^-^ (5) 

where G^*^ are the Sachs form factors and r^k denotes the third component of the isospin 
operator of the k-th nucleon. With our procedure the computational labour is reduced when 
the number of w-dependent form factors is reduced ^]. To this end we use the approximation 

GUQ') ^ G'^iQ'hiQl) (6) 

where 7(QL) = GUQl) / G'^iQl) , Ql = q' - and u^. = qV{2M). Similarly for the 
one-body spin current we use 

GW) ^ MQDGW) MQL) = (7) 



GUQ') ^ f^niQDGUQ') finiQl) = SwY- (8) 



For the usual dipole magnetic form factors, as used in this work, the above relations are 
fulfilled exactly and we have checked that the approximation provides a very good accuracy 
for G"^. In a future extension of our work to a high-g region, q > 500 MeV/c, we will 
use more sophisticated nucleon form factor fits. In these cases the above relations are only 
approximately fulfilled although we have checked that they still lead to excellent accuracy. 
The neutron electric form factor we use here is taken from [3] as used in [l^. With ([MH]) 
the one-body current is replaced by j*^-"^) G^^{Q'^)J^^^ where J^^-* is now given by 



J<"(q) 



k=l 



+ 



M 
i{ak X q) 



1 + r^k 
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The dominant contributions to the two-body current J*-^^ arise from the vr- and p-meson 
exchange currents. These currents are usually expressed in terms of "Seagull" and "true 
exchange" pieces. Thus we write here 



t(2) — ^ \p _i_ ; 
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ex ■ 



(10) 



We list in Appendix A the coordinate space representations of these currents with the 
corresponding values of coupling constants etc. Momentum space forms of these meson 
exchange currents are related to these coordinate space forms, apart from the multiplicative 
isovector electric form factor ^^(Q^) = (G^(Q^) - G^(Q2))/2, via 



ja(q)e 



iq-Rcn 



(11) 



where the super/sub-scripts above are those corresponding to the right hand side of ffTOj) . 
Finally we use the current operator J in the form 



(12) 



III. MULTIPOLE EXPANSION OF THE TRANSVERSE RESPONSE 

The dynamic calculations are performed in separate subspaces belonging to fixed angular 
momentum J and its projection M (see also [5]). One can account for M-dependencies 
analytically via performing a multipole expansion of Rt- To this end we use a decomposition 
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into multipoles of the transverse current. This decomposition shall also allow us employing 
an alternative expression for the transition operator, see below. The transverse current is 
represented as 

J^ = 4vr Yl E^'"^^-(^)Yf«l*(q)- (13) 

A=el,mag jm 

Here q = g^^q and are electric and magnetic vector spherical harmonics 



isj and e=0 



when A = el or e=l when A = mag. This then allows the transverse response to be written 

as 

RT{q,uj) = -^ J2 Y.^2J + l){Rryf (14) 

A=el,mag Jj 

where 

{RtYj' = ^ rf/(g^l,|M//(J,M))(vl/^(J,M)|g^l,)5(E^ - Eo - u), (15) 
J and M are the final state angular momentum and its projection, and l^'j^/) is given by 

\Ai) = ['^'^\MJo))]jM. (16) 

In Eq. ( |T5l) M is arbitrary. 



In terms of the more standard multipoles and vector spherical harmonics Yj^ we can 
write 

/ ■ . \ 1/2 / . N 1/2 

q-mag — rpj /-|Q\ 

jm jm \ / 

where 

TL = ^J d% (Y;„(q) ■ J(q, U)) . (19) 

Since charge has to be conserved it is well known that the above expression for TJ-^ can be 
rewritten as 

where pjm is a charge multipole of the charge density operator p defined by 

Pjmiq) = j dn^ Yjm{ci)p{(l). (21) 

We shall refer to the form of 7^-^ in f[T7l) as the direct form and to that in fl2U]) as the Siegert 
form. The first term of (1201) will be called Siegert operator, while the second term is the 
residual term. Appendix B gives the multipole operators Tj^ for the one-body currents 
while Appendix C lists them for the tt and p exchange currents. 
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IV. CALCULATION OF THE RESPONSE 



The techniques we use in calculating the response have been largely set out in [J]. Here 
we add some extra detail which arises in the case of the transverse response. The Lorentz 
transform of the partial response (i?^) of Eq. (1151) is given by 

The sum in (l22l) corresponds to transitions to discrete levels with excitation energy Un- In our 
A=3 case there exists only one discrete contribution corresponding to Ml elastic scattering. 
In (12^ the response is supplied with an additional superscript a. It specifies separate 
contributions to the response {Rt)j^ of Eq. flTBl) . e.g. a given a determines the isospin of 
the final state. In addition it specifies contributions that correspond to components of the 
multipole operators with different nucleon form factor dependencies. 

It was pointed out above that one-body and two-body currents have different u- 
dependence through their different form factors. Therefore we need to calculate the re- 
sponses with the individual parts of the current i.e. J^^^3^^\ J^-^^J^'^\ and J'^^^J*^^^ The 
corresponding partial response functions will carry the superscripts a = {11, 12, 22} so that 
the response (-Rt)'j' would be expressed as 

+A{GUQ'mRT{q,uj))f'\ (23) 

Additional cu-dependence of the electric multipole operators arises when they are used in 
the Siegert form, i.e. in the form of Eq. (1201) . Due to the additional cj-dependence of the 
first term in ( |20|) we calculate separately the response originating from this term (~ u"^), 
the response originating from the second term in fl20|) and the cross-term response (~ uj). 
For the same reason as in Eq. (!23|) each of these responses is in turn broken into the one- 
body piece, the two-body piece and the cross piece which are calculated separately. The 
superscript a in ( |22i) enumerates all these various cases, so that the response {RtY/^ is 
a sum of the responses [Rxyf'"' multiplied by products of nucleon form factors times u"', 
n = 0,l,2. 

As described in [3| the transforms are determined dynamically. In the present case the 
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transforms are obtained from 

The calculation fl24l) is M-independent and is performed in separate subspaces belonging to 
given isospin and parity. Parities are determined by the multipole order j and the choice of 
A=el/mag. For a given A, parity, and J only one value of j is possible in our case. 

To pass to responses one needs to invert the transforms. This may be done either sep- 
arately for each transform ^-^'^'^ using Eq. (!22l) or for their sums at the same a. One may 
define the responses i?^ = X]a=c1 mag -^t"' where (c.f. f|T^ ) 

4'"(g,o.) = -^5^(2J+ l)(i?T)f "(g,^). (25) 

One also defines the corresponding transforms 

$^'"(g, an. a,) = —— 5^(2J + l)$S'''^(g, a^, a,). (26) 

They are related to the responses (|25ll in the same way as in (!22l) . All the various are 
inverted separately to get -R^". Our inversion method and more information concerning the 
inversion can be found in jsl-ligl. I20I . We take a/ = 20 MeV and distinguish between the two 
isospin cases T=l/2 and 3/2, since the corresponding responses have different thresholds 
and thus the inversion can be carried out more precisely. After having inverted both cases 
we sum up the two results. For the magnetic part of the response we invert the Ml transition 
to the final state with J"^ = separately, since, as mentioned above, it contains an elastic 
contribution. This elastic contribution can easily be determined by choosing a very small 
value for aj and thereafter its effect on the transform can be subtracted leading to a LIT of 
a purely inelastic response. 

As mentioned earlier charge conservation leads to the equality of the Siegert and direct 
forms of the transverse electric multipole operator. This provides an important check on 
our procedures especially with respect to the implementation of the MFCs. The BonnRA 
potential contains more than just vr- and p- meson exchange but we expect that taking 
account of MFCs from only these two exchanged particles should lead to the dominant MFC 
contribution in our kinematical range, while additional MFC effects are partially taken care 
of by the Siegert operator. A good test for the implementation of the MFC is provided by 
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using a simple vr + p OBEP with their corresponding MECs. In this case charge conservation 
should be exact and the transverse response should be independent of whether one uses the 
Siegert or the direct form of Tj^^. We have made such tests at ^=10, 300, and 500 MeV/c 
and found very good agreement between the results of the two calculations 21 1. 



V. RESULTS AND DISCUSSION 



We have selected the momentum transfers q = 174, 250, 400, and 500 MeV/c for a 
calculation of Rxiq-iOj) in a large cu-range. In addition we consider the low-cj part of Rt at 
q = 174, 324, and 487 MeV/c for which cases we take a maximal value of J = 7/2. For 
the other g- values a different choice for J"^°-^ is made: 11/2 (g=250 MeV/c), 15/2 (g=400 
MeV/c), and 19/2 (g=500 MeV/c). We have checked that with these settings very good 
convergences of the multipole expansions of Rt are obtained in the requested energy ranges. 

In the discussion we compare results calculated with the various current operators of 
section II (both with direct and Siegert forms) representing the following contributions: (a) 
one-body, (b) one-body and implicit MEC via Siegert operator, (c) one-body, vr- and p- 
MEC, and (d) one-body, tt- and p-MEC plus additional MEC via Siegert operator. If exact 
charge conservation was satisfied then the results of the direct calculation (c) would agree 
with those of the Siegert form (d). 

In Figs. 1 and 2 we show the various current contributions to Rt- It is readily seen 
that there are rather strong MEC effects: 15-30 MeV above threshold MEC enhance Rt by 
more than 30% for the two higher g-values (very close to threshold even by up to 200%, 
see Fig. 4); they increase the quasi-elastic peak height by 10% (g=174, 250 MeV/c), 7% 
(g=400 MeV/c), and 6% (g=500 MeV/c); for lower q they also lead to large effects in the 
high-energy tail (e.g. at pion threshold: increases of 180% (g=174 MeV/c), 95% (g=250 
MeV/c), 22% (g=400 MeV/c), and 5% (g=500 MeV/c)). In general, relative contributions 
of MEC are determined mainly by distances \uj — Upeakl- This is natural since the peaks 
correspond to maximum contributions of one-body operators. 

It is also seen that Siegert contributions remain quite small in and below the quasi- 
elastic peak. On the other hand they become more important with increasing energy (e.g., 
at pion threshold and g=174 (250) MeV/c, enhancements are of 130% (55%) of the one- 
body contribution). In addition to the fact that in general MEC contributions are rather 
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small in the peak as compared to one-body contributions, Siegert contributions are strongly 
suppressed in and below the peak by the factor uj/q m f l20|) . The approximate transition 
operator we discuss takes account of MEC only via the Siegert operator, i.e. the charge 
operator from fl2U]) . As it is seen from Fig. 2, in the tail region this approximation provides 
the response rather close to the true one at the lowest q value q = 174 MeV/c. This agrees 
with the well-known fact that in moderate energy photodisintegration processes {u = q) 
MEC contributions are largely included by the Siegert operator. 

It is interesting to note that there exist additional Siegert MEC contributions beyond the 
TT- and p-MEC entering the direct calculation. This is due to the fact that vr and p exchanges 
constitute only the dominating part of a consistent exchange current with the BonnRA 
potential. Other two-body currents are induced by momentum and spin-orbit dependent 
potential terms. In addition also three-body currents, originating from the TM-3NF, could 
lead to Siegert contributions. Effects of the Siegert operator beyond vr- and p-MEC were also 
found in the proton-deuteron radiative capture with the BonnCD-|-A potential 221]. It is 
seen from Fig. 1 and Fig. 2 that such effects are small for energies far away from the photon 
point. Indeed, there the complete calculation via direct inclusion of MEC operators and the 
complete alternative calculation that involves the Siegert operator have led to results close 
to each other. However, closer to the photon point the additional Siegert contributions can 
lead to corrections of the order of 10%. 

In Fig. 3 we show our Rt results in comparison to experimental data. For q= 250 and 400 
MeV/c one finds good agreement. However, data are not precise enough to allow a definite 
conclusion about the MEC contribution. As opposed to the lower q cases we find at g=500 
MeV/c a difference between the theoretical and experimental peak positions. The shift 
amounts to about 5-10 MeV. Relativistic effects, in particular those arising from corrections 
to the kinetic energy, might be responsible for this difference. In fact in {g] it was shown 
for the longitudinal response function RL{q,uj) that such effects lead at g=500 MeV/c to a 
shift of the peak position by 6 MeV. 

In Fig. 4 we depict various Rt theoretical and experimental low-energy results at g=0.882, 
1.64, and 2.47 fm~^ corresponding to about 174, 324, and 487 MeV/c, respectively. We do 
not show the contribution of the Siegert operator, since, as shown in Fig. 1, its effect is very 
small at low energies. One sees that the MEC contribution can be very important, e.g. at 
q=487 MeV/c one finds an increase of about 200% close to threshold. Contrary to the cases 
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shown in Fig. 3 one can make a definite conclusion about the MEC contribution. It is evident 
that they lead to a considerably improved agreement between theory and experiment. For 
the two higher g-values theoretical and experimental results agree very well, whereas for 
g=174 MeV/c the theoretical result underestimates experimental data somewhat below 10 
MeV. A better theoretical description of the g=174 MeV/c data is found in [l], where the 
AV18 NN potential and the UrbanalX 3N-force is used as nuclear interaction. However, the 
Coulomb force was not included in the final state interaction. The effect of such a neglect 
is illustrated in Fig. 4 for the case in discussion. Within 2 MeV above threshold it leads to 
an increase of more than 10%, while at 5 MeV above threshold the effect still amounts to 
4%. In this way the theoretical results are shifted closer to the experimental data, but the 
effect is too small to reach a good agreement at low energies. 

We summarize our results as follows. We have calculated the transverse form factor 
RT{q,uj) considering besides one- and two-body currents also the so-called Siegert operator. 
As nuclear interaction we have taken the BonnRA NN potential and the Tucson-Melbourne 
TM' 3N-force. Since we are particularly interested in the MEC effects and the role of the 
Siegert operator we have chosen the BonnRA potential, for which the important vr- and 
p-exchange currents are directly determined by the potential model. It is true that also 
for more phenomenological NN potentials, e.g. AV18, a consistent vr- and p-MEC can be 



constructed 



la, [iTI], but to this end one has to interpret the isovector part of the 



phenomenological potential as an effective vr- and p-exchange. 

We find that MEC provide very strong contributions both at lower energies and in the 
high-energy tail while giving a moderate increase to the height of the quasi-elastic peak. 
Siegert contributions are unimportant in and below the quasi-elastic peak. They become 
considerably more sizeable at higher energies, but additional MEC contributions have also 
to be taken into account and thus a calculation, where in addition to the one-body current, 
MEC currents are taken into account via only the Siegert operator is not sufficient. On the 
other hand also a calculation with only one-body currents plus vr- and p-MEC may not be 
sufficient at higher energies, since, as we have shown, effects due to additional two- and three- 
body currents can become important. To include at least a part of these additional exchange 
effects it is better to work also in this case with the Siegert operator. The appropriate place 
to study the structure of MEC is the energy region below the quasielastic peak. Indeed, 
the contributions of MEC are large in this region while they are rather small in the peak. 
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and beyond the peak they are partly represented by the Siegert operator, i.e. the charge 
operator. 

In comparison to experiment relatively good agreement is obtained at g=250 and 400 
MeV/c, while at g=500 MeV/c, presumably because of relativistic effects, the position of 
the theoretical quasi-elastic peak is located somewhat above the experimental one. Close to 
threshold one finds very strong MEC contributions. They are necessary in order to achieve a 
good description of the experimental data at q=324 and 487 MeV/c. Also at g=174 MeV/c 
they lead to an improved agreement with experiment, but in the range from threshold to 5 
MeV above the theoretical result underestimates data somewhat. 

In future we plan to investigate the momentum range 500 MeV/ c < g < 1 GeV/ c consider- 
ing relativistic corrections for the one-body current operator and performing the calculation 
in a reference frame where relativistic effects in the kinetic energy are minimized |6|. We 
also plan to study isobar current contributions including A(1232) degrees of freedom. 
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APPENDIX A: CONFIGURATION SPACE vr AND p MECS 



For convenience we list below the well known vr and p configuration space exchange 
currents. 



[{a, ■ V0a,-5(x - r,) - (a, ■ V,)a.5(x - r,)] J] |r, - r,|), (Al) 

k=l 

^ hi [Y(^l, |x - r,\)VJ'^^Ll, |x - r.l) - Y(^l, |x - r,\)V^Y(^l, |x - r,|)], (A2) 

3 

[{ai X V,) X a,-5(x - r,) - (a, x V^) x ^^^(x - r^)] K^i^l - rj)!)- (A3) 



3 

X 

fc=l 

2 



A:=l 



3 

X 



K Ix - r,|)V.r(/x^, |x - r,|) - |x - r,|)V.y(/.^, |x - r, |)] . 



k=l 



(A4) 



Here Y{m,r) = e "^"^ jr. We list the coupling constants, the masses /x^, and the regular 



ization constants /i^, where a = vr or p, taken from ll|: 



/o = ^Inn = 0.0805, i = 1.2, ^ = 6.1, 
4vr 47r 

/i? = m^, p^ = A^ + 10 MeV, p.'^ = Aa - 10 MeV, (A5) 



A^ = 1.3 GeV, Ap = 1.2 GeV, (A6) 
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APPENDIX B: MULTIPOLES OF ONE-BODY CURRENTS 

In the following the non-relativistic expressions of the electric and magnetic multipoles 
for the one-body currents are written. Each of them is decomposed in a convection and a 
spin current. For the magnetic multipoles one has 



with: 



|-^Ji-.(4'-D|ij-i(i-;) ® ".U \ . (B2) 



rp3, corw _ _;_ 



1 /1 + 7 1-7 



+ 



(B3) 



The quantity 9^ is defined by the relationship —id'^ = p'^. If the last Jacobi vector is defined 
as U-i = ViA-l)/A [ta -{A- 1)-^ ^fj,^ r,] then 



A-1 
A 



1/2 



d 



A-l,ti 



Similarly we write the one-body multipoles contributing to 7^'^ as 



rpl \ ^Irpl,Spm/ ■\ . rpl,COnv / 



(B4) 



where / = j ± 1. One obtains 



T; 



ijitljSpin, 



' M2 



and 



T; 



ijitl, conv , 



t) = ±- 



+ 



1 /1+7 1-7 



' ^2jll^^ ^^(9^'^ ® (B5) 



(B6) 



The term proportional to jj{qr'j) in (1B6P above cancels when one forms the electric multipole 



16 



APPENDIX C: rfj^ MULTIPOLES OF vr AND p MECS 



Here the Tjj^ multipoles are given for the "12" pair. The total result should be multiplied 
by 3 to account for three pairs of identical particles in the trinucleons. 

1. TT- Seagull 

^ ( ' 

Ipa a'C 

1 a \ f a' £ L \ \ L i a' \ \ L 1 C 
000l\000l\aCl\\pJl 



T 



JM 



(—)' '''''' ^-^y^^^^'l^ + {-^Y^'Y^pLaa't 



(n X T2), . 



(Ci) 



Here r = r2 — ri, z = — [(ri + r2)/2 — Rem], f = r/r, z = z/z, {H^{r))' = dH'^ {r) / dr , and 
H'^{r) = X]fc=i ^fc ^(/^L ''^)) where the constants are given by Eqs. (lASp . (]A7|) . We denote 
the spin-coupling [cti C?) 0-2] p^m by T,'^^. 

2. 7r-Exchange Current 

The multipole for the true-vr exchange is the sum 

JM ~ -'■ JM ' JM + 7 



■ JM 



- JM 



■ JM 



where 



-lT{lfpL {L'faa'C 




"12 



M 



in X rs). , (C2) 



pL,X2 
■ JM 



4vr q' f fo 




epaa' L'C 

1 £ a\ ( L' p L 

oooyl^ooo 

L J M 




(C3) 
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T 



47r4v^g / /o 



JM 



1 i J' \ 1 J' a \ L 1 L' \ L' i a' \ \ L' £ a' \ \ 1 £ J' 
OOOyl^OO oJyOOoJyOOoJ]^aCf\\alf 

L' f 
1 1 1 



}j.'{q^)^^'kq,r) [n.(z)®n(f)]^®sw (nxTs),. (C4) 

M 



L 1 J 



Here r) = ELi C;(? 

,r,jj,l), where the functions (j)^i{q,r,m) are defined in [27|. 



^3 7 TT 



The multipoles ( 1C1I - IC4I) are reaL 
3. p- Exchange Currents 

The multipoles of the p-exchange currents can be obtained from the above 7r-exchange 
currents by means of the following replacements: /i^ ^ /i^, hi ^ h'j^, 

2 1 n / r \ 2 

(C5) 



mj 47r \2MJ \ gj 



and by inserting in each equation above the factor 



-6(-ir 



1 1 1 
1 p 1 



(C6) 
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CAPTIONS TO FIGURES 



FIG. 1: Effects of the various contributions on Rt in the quasi-elastic region at g=174 (upper left), 
250 (upper right), 400 (lower left), and 500 MeV/c (lower right): one-body (dotted), one-body 
-|- implicit MFC via Siegert operator (dashed), one-body -|- vr-MFC -|- p-MFC (dashed dotted), 
one-body -|- vr-MFC -|- p-MEC -|- additional MFC via Siegert operator (solid). 



FIG. 2: As Fig.l but for the high-energy region. 



FIG. 3: Comparison of theoretical and experimental Rt at g=250 (upper panel), 400 (middle 
panel), and 500 MeV/c (lower panel). Theoretical Rt with contributions: one-body (dotted) and 
one-body -|- vr-MFC -|- p-MFC -|- additional MFC via Siegert operator (solid). Experimental data 



from ^] (triangles), [24| (circles), and 



25l | (squares). 
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FIG. 4: Comparison of theoretical and experimental low-energy Rt at g=0.882 (upper panel), 1.64 
(middle panel), and 2.47 fm~^ (lower panel). Notation of curves as in Fig. 3, but additional curve 
in upper panel for total result in case that Coulomb force is neglected in the final state interaction 
(dash-dotted). Experimental data from 26l |. 
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